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SUMMARY 


A nonlinear theory for sound propagation in variable area 
ducts carrying a nearly sonic flow is presented. Linear acoustic 
theory is shown to be singular and the detailed nature of the 
singularity is used to develop the correct nonlinear theory. 

The theory is based on a quasi-one dimensional model. It is 
derived by the method of matched asymptotic expansions. 

In a nearly choked flow the theory indicates the following 
processes to be acting: A transonic trapping of upstream 

propagating sound causing an intensification of this sound in 
the throat region of the duct; generation of superharmonics and 
an acoustic streaming effect; development of shocks in the 
acoustic quantities near the throat. 

Several specific problems are solved analytically and 
numerical parameter studies are carried out. Results indicate 
that appreciable acoustic power is shifted to higher harmonics 
as shocked conditions are approached. The effect of the throat 
Mach number on the attenuation of upstream propagating sound 
excited by a fixed source is also determined. 


I. INTRODUCTION 


The acoustic behavior of variable area ducts has recently 
been the subject of much attention in connection with sonic 
engine inlets. Observations of a correlation between axial Mach 
mxtnber and attenuation of sound radiated upstream from such 
inlets has led to a realization that nearly choked inlets can 
be effective in suppressing upstream propagating sound Ref. 

1-3. However until recently. Ref. 4,5 , there has been no 
parallel theoretical study to explain the physical processes 
that are responsible for this attenuation. In Refs. 4,5 a 
new nonlinear theory describing sound propagation in ducts 
having a near sonic throat section has been derived and some 
numerical results given. The present report will give a detailed 
derivation of the nonlinear theory and present analytical and 
numerical studies to illustrate the predictions of the theory. 

The theory provides a physical and quantitative understanding 
of the propagation process in the throat region. The nonlinear 
interaction between the sound field and the basic flow in the duct 
is shown to give rise to: 

(1) An intensification of the upstream propagating sound 
in the throat region due to a transonic type trapping; 

(2) The generation of higher harmonics of the fundamental 
(source) frequency; 

(3) An acoustic streaming effect; 

(4) The development of shock waves in the perturbation 
quantities in general. 

A condition on the source strength and frequency/ ^nd the 
throat Mach number which determines whether a shock will occur 
has been derived and numerical results have been obtained which 
illustrate the development of such "acoustic" shocks. All 
numerical results presented in this report are for shock-free 
solutions. These results show that a large fraction of the 
acoustic power transmitted through the duct can be shifted into 
the superharmonics even in the absence of shock formation. The 
attenuation for upstream transmitted acoustic power as a function 
of the throat Mach number for a fixed source strength is also 
studied. The nonlinear theory predicts that the acoustic power 
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should be of the order of e where e is the deviation of the 
throat Mach number from unity while a naive linear theory predicts 

an dependence. 

The nonlinear theory is based on a quasi-one dimensional 
model which yields results of practical interest and provides 
basic understanding of the physical phenomena involved while 
remaining relatively simple to treat analytically and numerically. 
The theory is derived by observing that linearized acoustic 
theory fails to descirbe propagation of sound in ducts containing 
a near sonic throat section; the failure being manifested in 
singular behavior (infinite amplitudes) of the acoustic quantities 
in the throat region. A discussion of the failure of linear 
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theory is presented in Section IIB. A more complete treatment 
is given in Ref. 6, In Section IIIA the localized nature of 
the singularity is exploited in deriving the correct nonlinear 
acoustic theory by use of the method of matched asymptotic 
expansions. The theory is applied to yield numerical results 
for several cases of sound propagation from an acoustic source, 
located at and upstream or downstream of a near sonic throat. 
These results are discussed in Section IIIE. 

The qualitative basis of the theory is relatively easily 
understood in terms of the Riemann invariants of the acoustic 
field. If P(x,t) and Q(x,t) represent the downstream and 
upstream propagating acoustic waves, respectively, then the 
results of the theory are as indicated in Figure 1. Nonlinear 
transonic effects occur in the lowest order acoustic perturba- 
tion quantities whenever the perturbation level in the throat 
region is of the order of magnitude of the deviation of the 
throat Mach number from unity. This can occur in two 
ways: If the acoustic source is at the throat and has a 

strength of 0(e) , or if the acoustic source is outside the 

2 

throat and has a strength of 0(e) where e is the Mach niimber 
deviation. To fix ideas consider the latter case, illustrated 
in Figure 1. Directly upstream of the source a continuous 

reflection process occurs because of the variable area, so 
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that both P and Q waves are present at 0(e ) and are described 
by linear theory. However, the asymptotic inner solution shows 
that the Q wave is strengthened to 0(e) by nonlinear interactions 
with the steady flow which occur when x = 0(e) . Since the inner 
region is small compared to the scale of variations of the steady 
flow, no reflection of the Q wave occurs at 0(e) : The P wave 
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remains at 0(e ), the relationship between them being ultimately 
determined by an acoustic boundary condition applied at some 
X < 0. Thus, as might be anticipated on physical grounds, the 
downstream wave can be considered to propagate through the throat 
essentially unaffected by the near-sonic nonlinearity, while the 
upstream wave is intensified in the throat region due to a trans- 
sonic trapping effect. The trapping of the upstream wave occurs 
because its velocity relative to a fixed observer is practically 
zero. It thus has a long residence time in the throat region 
of the duct. 

Other recent work on acoustic propagation in high subsonic 
Mach number flow has been numerical, Ref . 7,8, In both cases 

a quasi-one dimensional model is treated. In Ref. 7 an 
expansion in harmonics has been used and the resulting coupled 
equations integrated numerically. Some results are presented 
only for unshocked solutions and for a special reflection 
coefficient. They support the theoretical conclusions of the 
present paper. In Ref. 8 the initial value problem is inte- 
grated to a steady state and some numerical results are presented 
for the acoustic velocity potential. There is no discussion of 
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the harmonic content of the solutions and the question of 
shocks is not considered. 

Many earlier authors have studied linear quasi-one dimen- 
sional duct acoustics, but, in general, these studies have 
been concerned neither with the behavior nor with the validity 
of the linearized theory as the mean flow approaches sonic 
speed. A comprehensive review of earlier work can be found 
in Ref. 9, Tsien, and later Crocco Ref. 10, derived an 
exact solution which is discussed in Appendix D of the present 
report. Eisenberg and Kao Ref. 11 found another exact 
solution corresponding to a special, but practically unreal- 
istic, area variation. Numerical evaluation of their analytical 
results clearly exhibits singular behavior as the axial Mach 
number approaches unity. Davis and Johnson Ref. 12 performed 
numerical computations of the linearized solution for certain 
acoustic boundary conditions, while King and Karamcheti Ref. 13 
employed the method of characteristics to obtain similar numeri- 
cal results. 

The earlier numerical studies cited above yield results 
for sound propagation through a throat so long as the Mach 
number there is sufficiently far from unity. They include no 
attempt to resolve the specific behavior of the linear solution 
in the event that the throat flow approaches sonic speed, 
although they are certainly useful in indicating that difficul- 
ties exist with the theory in such a circumstance. Several 
numerical solutions illustrating the development of the acoustic 
singularity were presented in Ref. 6. 

The research presented in this report is the result of a 
joint effort by the author and Professor M. K. Myers of the 
George Washington University, Joint Institute for Advancement of 
Flight Sciences. Professor Myers' work was supported by NASA 
Langley Research Center under NASA Grant NGR-09-010-064 . 
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II. SYMBOLS AND NOMENCLATURE 


a parameter measuring throat curvature 


A(x) 

area of 

duct, m^ 

+ 

+ 



AT , 
In' 

B7 

In 

see 

Eq. (2.41) 

+ 



+ 

A~ 

n 

magnitude 

of At 

In 

A*, 

B*, 


see Eq. 

B 

see Eq. 

(2. 

1) 


B(x), Bj^ see Eq. (2.34) 

Cg stagnation sound speed, ms~^ 

c (x) steady flow sound speed 

see Eg. (2.33) 

+ + + 

^02' *^11' ^^12 (2.42) 

+ + 

(x) , (x) inner Fourier coefficients 

d, e geometric parameters in area variation 

®ln' ^2n 

it + + 

E^, real and imaginary parts of I^ 

Ff, F 2 , F^, F^ Hypergeometric functions 

g(x;e) see Eq. (3.29) 

G(x;e) see Eq. (2.5) 

Gq G (x; e = 0) 

H(x) see Eq. (2.30) 

+ 

I^ see Eq. (3.32) 

j, k, £ parameters in Crocco-Tsien solution 
J see Eq. (3.19) 

K parameter in Crocco-Tsien Mach number distribution 



L characteristic length associated with area variation^ m 
M(x;e) steady flow Mach nximber 

Mq ( x) M(x;0); leading order Mach number 

M^ (x) , M 2 (x) see Eq. (2.12) 

^10' ^20 local expansion Mach number terms 

m^ (X) , m^ (X) inner expansion Mach number terms 
Nq, N source amplitudes 

P(x,t) downstream perturbation Riemann invariant 

P^, pi first and second order downstream perturbation 
Riemann invariants (inner region) 

_ -2 
P(x,t) total fluid pressure, kgm 

Q(x,t) upstream perturbation Riemann invariant 

oi, oi first and second order upstream perturbation 
Riemann invariants (inner region) 

r(x,t,e) outer, linearized density 

r_(x,t), r, (x,t) first and second order terms of outer 

density 

+ 

r^^(x) complex amplitude of density 

r^, r^ first and second order terms of inner density 
R^ stagnation density 

s see Eq. (2.11) 

Sq acoustic source strength 

t time, second 

T see Eqs. (2.26) 

u(x,t) total fluid velocity, ms ^ 

U(x) basic steady flow, velocity, ms ^ 

u(x,t) perturbation velocity 

Ug / U source amplitudes, ms ^ 
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V 


see Eqs. (2.26) 

W see Eq. (2.29) 

W“, W^, W* see Eq. (2.41) 

X axial distance along duct 

X stretched distance (inner variable) 

x^ boundary location 

z Crocco-Tsien axial distance variable 

a see Eq. (A7) 

a (x) duct area ratio 

6 dimensionless Crocco-Tsien frequency 

Y ratio of specific heats 

see Eq. (D.IO) 

r see Eq. (A7) 

6 outer acoustic disturbance level 

A Crocco-Tsien parameter 

^ throat Mach number deviation from unity 

0 “ phase angle inner solution 

n 

0 , 0 phase angles outer solution 

rn yn ^ ^ 

A, Aj^, A^ eigenvalue of linearized acoustics local solution 

y(x,t;e) outer linearized perturbation velocity 
yg(x,t), y^(x,t) see Eq. (2.22) 

y^^(x) complex velocity amplitude of n-th harmonic 


7 



see Eqs. (3.9) 


y^(X,t), y2^X,t) 

V see Eq. (3.32) 

IT see Eq. (3.15) 

p(x,t) total fluid density 

p(x,t) perturbation density 

T characteristic parameter inner region 

(j)(t) throat source function 

(jo radian source frequency (radians per sec.) 

dimensionless frequency 


Special Notation 


indicates absolute value of a real quantity or 
magnitude of a complex quantity 


subscript n 


indicates a 

superscript 

+ 

indicates 

superscript 

i 

indicates 

superscript 

0 

indicates 


harmonic 

X > 0, X < 0, respectively 
inner variable 
outer variable 
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III. FORMULATION AND OUTER EXPAtJSIONS 
A. Perturbation Equations 


Consider the propagation of sound in a variable area 
duct carrying a homentropic inviscid ideal gas flow. The 
acoustic wave lenqth is assumed sufficiently large and the 
area variation sufficiently slow, that the field can be 
described by the equations of quasi-one dimensional gas 
dynamics (Ref. 14) 

°s^t ° 

CgU^ + uu^ + Px = 0 ' (2.1) 

P 

= constant = B 


In equations (2.1) p, p, and u are the total fluid pressure, 
density, and axial velocity, and A(x) is the duct cross 
sectional area. The dimensionless independent variables x 
and t are measured in units of L and L/c^ respectively, where 

L is a characteristic length associated with the area variation 
and c is the stagnation value of sound speed in the gas. The 
geometry of the problem is as indicated in Figure 2 where the 
origin of x corresponds to a throat; A' (0) = 0. 

If the velocity and density in the basic steady flow in 
the duct are denoted by U(x) and R(x) respectively, then from 
( 2 . 1 ) , 


UR' + RU' + RU(^) = 0 
UU' + yBR^~^R’ = 0 , 


( 2 . 2 ) 


where the energy relation in Eqs. (2.1) has been used to 
eliminate the pressure from the system. In order to study 
solutions to the system (2.1) which are small perturbations 
about the steady values U and R, it is convenient at the outset 
to define dimensionless variables u(x,t) and p(x,t) by 

u(x,t) = U (x) [l+u(x ,t) ] , p(x,t) = R(x) [1+p (x, t) ] , (2.3) 


where it is assumed that U(x) 0. It should be noted that 
the usual acoustic velocity and density are given by U(x)u(x,t) 
and R(x)p(x,t), respectively. 



Substituting Eq. (2.3) into Eq. (2.1) and employing the 
steady relations (2.2) yields the system of equations on u 
and p in the form 

+ M[(l+u)p^ + d+p)u^] = 0 

G^/^u^ + M(l+u)u^ + ^ (l+p)^"^p^ 

+ ^ [ (1+u)^ - (1+p)^"^] = 0 

In Eqs. (2.4), M(x) is the flow Mach number U(x)/c(x), c (x) 

2 Y-1 

is the speed of sound in the steady flow (c = yBR ) , and 

G(x) = (^) = 1 + , (2.5) 
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(2.4) 


the latter expression following from the Bernoulli relation 
implied by the second of Eqs. (2.2). 

For the purposes of physical interpretation and ease of 
asymptotic matching, the system (2.4) is recast in Riemann 
invariant form (Ref. 14) . The upstream and downstream pro- 
pagating perturbation Riemann invar iats are defined by 


P ( X , t ) = Mu + 


Q(x,t) = Mu - 


Y 


Y 


^0 


(1 + p) 


(1 + p) 


Y-I 

2 


Y-1 

2 




( 2 . 6 ) 


Y-3 


If the first equation in system (2.4) is multiplied by (1+p) 
and the second by M(x), then adding and subtracting the 
resulting two equations and using Eq. (2.6) and analogous 
expressions for the partial derivatives of P and Q yields 




Y-1 


M(l + u) + (1 + p) 


Y-1 


J^uM(l + u) + (l + p)^ =0 , (2.7) 


M 

a/2 
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Y-1 

|M(l + u) - (1+ p) ^ [(1 + u)^ 

Y-1 


;x72 [» 


-j^uM(l + u) - (1 + p) ^ 


To complete the elimination of u and p the relations 


(l + U) ‘ 1 . 


Y-1 




(1 + P) ^ ^ 

are substituted in 

Eqs. (2 

^t + 

1 

G^/2 

M + 

1 p 

+ 

MM' 

g^/2 1 

[P+Q 

_M 

, (P+Q)^ 
2 

4M 


M' 

fP+Ql 

M + 1 + 


G^/2 

2M 

^ y 

Qt + 

1 

M — 

1 + 3-Y . 

g"^2 

i*l 

J- + 4 

1 

MM' 

■p+Q 

. (P+Q)^ 

+ 


_M 

4m2 


M' 1 

r ^ 

P+Q 

[m - 1 + 


?T2 i 

2M 

^ y 


1 + (P - Q) 


P + ^ Q|P, 




(p-o) - (p- 




^ = 0 


P + ^ QIQ, 


=] 


(P-Q) - (p. 


P + Q 
4 4 


^ = 0 


(1 + p)^"^] 

(2.7) 


0,^] 

( 2 . 8 ) 

0)^] 
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It is assumed that the basic steady flow quantities U 
and R are known from solving the system (2.2) subject to 
suitable boundary conditions applied far upstream and down- 
stream of the throat of the duct, x = 0. Then Eqs. (2.4) or 
equivalently (2.8) are a quasi-linear system of partial 
differential equations with strong spatial dependence in the 
coefficients (due to the steady flow Mach number terms) . 

These equations describe perturbation solutions of any size 
about the known steady flow. 

Since propagation of sound in near sonic flows is the 
problem of interest, it will be further assiamed that the 
boundary conditions on the basic flow yield a throat Mach 
number which is near unity- Hence, a parameter e = l-M(O) 
is introduced into the analysis, and the Mach number M is 
regarded as a function of e as well as of x. In the analysis 
to follow, e will be assumed to be positive, i.e., the basic 
flow is subsonic throughout the duct. The coefficients in 
the system (2.4) or (2.8) then all become functions of e, and 
the nonlinear theory will result from considering these per- 
turbation equations for small e. In order to analyze these 
equations the detailed dependence of the basic steady flow 
on the parameter e must be found. This is discussed in the 
next section. 


B. Basic Steady Flow 

The elementary equations of quasi-one dimensional steady 
flow (2.2) are discussed in detail in numerous texts; for 
example, a particulary comprehensive treatment has been given 
by Crocco (Ref. 15) . It is straightforward to express any of 
the fluid quantities in terms of the duct area A(x) or, 
equivalently, in terms of the Mach number M(x;e). This was 
first analyzed in Ref. 6. It is the purpose of the present 
section to determine the behavior of M explicitly as a function 
of X and e. 

To begin consider the well-known relation implied by 
Eg. (2.2), 


M' = 


-MGA' 


(2.9) 


(1 - M )A 

which becomes, after integration, 

a^M^(x)[l + ~ - (P -I] = M^(0)[1 + (2.10) 


where in Eq. (2.10) the definitions 
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and 


( 2 , 11 ) 


A(x) 

A(0) 


a (x) 


/ 


_ 2(y-1) 
(Y + 1) 


have been used. Equation (2.10) is an implicit solution of the 
differential equation (2.9) and hence an algebraic equation 
for M(x) . Figure 3 shows sketches of typical integral curves 
of Eq. (2.9) for the type of area variation to be used here. 
Curves such as AB in the figure, for which M remains less than 
unity for all x are of interest. Since l-M(O) = e is assumed 
small, it is natural to seek an expansion of M(x;e) in the 
form 

M(x;e) = Mq (x) + (x) + e^M 2 (x) + ... , (2.12) 

where Mq( 0) = 1. Substituting Eq. (2.12) into Eq. (2.10) and 
equating like powers of e one finds that Mq( 0) must satisfy 


(Y + 1) ct^M® 


= 1 + 


(Y - i)m; 


while (x) = 0, and 


M 2 (x) 


-2Mq[1 + 


(Y - DM 


(Y + 1) (1 - Mq) 


(2.13) 


(2.14) 


Obviously the expansion (2.12) is not uniformly valid near 
x = 0: The third term is of order e whenever 1 - M^ (x) is as 

small as e, or equivalently M 2 (x) becomes infinite as x ->■ 0. 

It remains to find Mq ( x) = M(x;0) in terms of x: i.e., 

to solve Eq. (2.13). It is assiamed for simplicity that a(x) 
is an analytic function near x = 0 with a' (0) = 0, although 
the analysis to be carried out is easily extended to cover 
cases where a may have jump discontinuities in its higher 
derivatives at x = 0. Thus one may write a(x) in the form 


a (x) 


1 + ax + bx + dx + ex + 


(2.15) 
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The nonlinear acoustic theory to be studied below is derived 
for ducts with area variations satisfying the condition 
a = A"(0)/2A(0) > 0. This is not a significant restriction 
since if a < 0 the duct cannot have a throat and the case 
a = 0 is obviously special for an actual duct. Expanding 
Mq (x) in a power series of the form 


2 

Mq (x) == 1 + + ... / (2.16) 

substituting Eqs. (2.15) and (2.16) into Eq. (2.13) and 
equating to zero corresponding powers of x, yields for the 
branch (x) £ 1, 

1/2 

Mq(x) = 1 - [ ^ 1^1 + "^20^^ • (2.17) 


Hence, the leading term of M(x;e) behaves as a piecewise 
linear function of x near the throat so long as a / 0. 
The actual value of m 2 Q in Eq. (2.17) is 


m 


20 


— + 


^ 2 ^ 2 


(5y - 3) a 
12 


the ± sign holding for x J 0 respectively. All the basic 
steady flow quantities can be expressed in terms of the Mach 
number 


M(x;e) = 1 - 


1/2 


"^ 20 ^ 


+ e M2 (x) + 


(2.18) 


where etc., depend on Mq (x) and a(x). This expansion 

is to be considered asymptotically valid as e ^ 0 for fixed x. 

It is nonuniforra near x = 0; the appropriate Mach number dis- 
tribution near x = 0 will be given in Section III— A. 


C. Singularity in Linear Acoustics 

In this section linearized acoustic theory shall be 
analyzed. In particular the breakdown of the theory for pro- 
pagation in high subsonic flow will be demonstrated and the 
detailed nature of the singularity in acoustic quantities 
shall be determined. These results are essential for the 
derivation of the correct nonlinear theory. 

If it is assumed that the perturbation quantities u and 
p are small then the system (2.4) can be linearized. Let 6 
be a small parameter measuring the order of magnitude of the 
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perturbations. The physical interpretation of 6 in terms of 
an acoustic distrubance level will be discussed in Section 
III-C. Expanding u and p in powers of 6 as 


u(x,t;e,6) = 6y(x,t;e) + o(5 ) + 


p(x,t;e,6) = 6r(x,t;e) + o(5 ) + ... 


(2.19) 


substituting Eqs. (2.19) into Eqs. (2.4) and neglecting all 
but first order terms in 6 yields the linearized acoustic 
equations 


+ M(r + y ) = 0 

t X X 

G^'^^y, + My + ^ r + [2y -• (Y~l)r] = 0 

^t X M X Q 1. V I / J 


( 2 . 20 ) 


Similarly, utilizing Eqs. (2.19) in Eqs. (2.6), the Riemann 
invariants are written in terms of the first-order pertur- 
bation quantities as 


P(x,t;e,6) 


5 (My + r) 


+ 0 ( 6 ^) 


Q (x, t; e, 6) 


6 (My - r) 


+ 0 (<S^) 


( 2 . 21 ) 


The linearized version of the system (2.8) is simple to 
derive, but it will not be necessary to use it in the analysis 
which follows. 

Equations (2.20), subject to appropriate boundary 
conditions, generally must be solved numerically because of 
their variable coefficients. It is necessary for present 
purposes to analyze the behavior of solutions to Eqs. (2.20) 
in the vicinity of the throat of the duct when the throat Mach 
number M(0) is close to unity. It is well known that the 
system (2.20) is singular at any point where M(x) = 1. This 
can be seen most simply by subtracting the two equations; the 
resulting equation has no y term, and the coefficient r 

becomes (M'^ - 1)/M, which vanishes as M ->■ 1. This can only 
occur at x = 0 for the duct of Figure 2 . The singularity at 
X = 0 implies that, in general, the acoustic quantities r and 
y will be singular when the flow is sonic there. Thus, as 
will be seen in what follows, r and y generally become 
arbitrarily large near x = 0 as M(0) approaches unity, there- 
by violating the assumptions made in deriving Eqs. (2.20) 
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that ]j, y , r and r all remain bounded. 

As a result of the singular behavior of the system (2.20) 
for high subsonic Mach numbers in the throat region, linearized 
acoustic theory fails for any 6, no matter how small, if the 
throat Mach number is sufficiently close to unity, i.e., for 
e small enough. In order to describe sound propagation in the 
duct in this event, one must reformulate the perturbation 
scheme to take into account non-linear terms in the system 
(2.4) which were neglected in Egs. (2.20). However, in order 
to proceed in this direction it is necessary to know precisely 
the nature of the singular behavior of the solutions y and r 
to Eqs. (2.20). This behavior has been recognized, but never 
resolved, in previous treatments of the system (2.20) Ref. 9. 

Since the steady flow depends on the parameter e = l-M(O), 
the coefficients in the acoustic equations and hence the 
acoustic quantities y and r are functions of e. For e << 1 
we look for solutions of the acoustic equations in the form 


r = rQ(x,t) + er^(x,t) + ... 
y = yQ(x,t) + ey^(x,t) + ... 


( 2 . 22 ) 


Inserting Eq. (2.22) into Eq. (2.20) and using expansion 
(2.12) for the coefficients, one obtains, after neglecting 
higher order terms in e 

t XX 

(2.23) 

M' 

^0^ ’^O. ''' ^O’^O mT ^ ^^^0 " ° 

t X U X u 

where is defined by Eq. (2.5) with M = M^ . M^ (x) is the 

leading order term corresponding to any Mach number distri- 
bution which yields a sonic velocity at the throat and a 
subsonic velocity throughout the remainder of the duct. These 
equations are singular at x = 0 since Mq( 0) = 1. Analytical 

solutions of the system for arbitrary M- (x) and arbitrary time 
dependence are not readily found. However, for harmonic time 
dependence the system can be reduced to a system of ordinary 
differential equations with a singular point at x = 0 and no 
other singular points within the duct. The nature of this 
singular point determines the singular behavior in the time 
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harmonic acoustic quantities. Explicit analytical results 
concerning the exact nature of the singular point and the 
dependence on (x) are now derived by use of series solution 

methods for linear ordinary differential equations. 

Thus solutions of Eqs. (2.23) in the form 

, /V inf^t 

r (x,t) = r (x)e 

^ (2.24) 

, . , , . inQt 

yQ(x,t) = Ug^(x)e 

are sought, where fi = o)L/c is the dimensionless fundamental 

s 

frequency and co is the radian frequency (fundamental) of the 
harmonic solution. The parameter n is introduced since super- 
harmonics of the fundamental will be shown to occur due to 
nonlinear effects. Substitution of Eq. (2.24) into Eqs. (2.23) 
and rewriting the resulting system in standard form yields 


M, 


On 


1 - M 


2 


M 


0 


On 


1 - M 


2 On '^On 

0 


(2.25) 


where 


S(x) 

T(x) 

V(x) 





2M* 


(Y - 


(Y - 1)M^ 


(2.26) 


As has been seen in Section II-B the behavior of 
1 - Mq ( x) near x = 0, and hence the nature of the singular 

point at X = 0, depends upon the corresponding behavior near 
x = 0 of the area A(x). The method of Frobenius (Ref. 16) 
will now be used to analyze the behavior of the general solution 
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to Eqs. (2.25) near x = 0 for any area variation with A" (0) 

^ 0. The case of A"(0) = 0 is discussed in detail in Ref. 6. 

Since A"(0) 0 the corresponding Mach number distribu- 

tion is given locally by Eg. (2.18). The Mach number is a 
piecewise linear function of x; thus one can analyze the 
system separately for x < 0 and x > 0. General solutions 
valid for x < 0 and x > 0 can be constructed by analytically 
extending the coefficients of the system to x > 0 and x < 0, 
respectively. This is done by extending M_ (x) as follows: 

To extend to x > 0 take 


Mq(x) = 1 + [iX±lI^] + 


and to extend to x < 0 take 


Mq(x) = 1 - [ + 


(2.27) 


(2.28) 


Consider the system (2.25) with M^ (x) given by either of 

Eqs. (2.27) or (2.28). In each case, S, T, and V are 
extended across x = 0 as analytic functions, and the only 
singularity is easily seen to be a simple pole by using 
Eg. (2.27) or (2.28). Hence both of the extended systems 
have a singular point of the first kind at x = 0 , and 
solutions can be found by the method of Frobenius. 

One accordingly rewrites the system (2.25), using matrix 
notation, as 


W = 


^0 
1 - M 


0 


” T -S" 

w , 

w = 

^On 

_-V S _ 



_ ^0n_ 


(2.29) 


As is customary in dealing with such singular points, one 
exhibits explicitly the simple pole by writing 


1 - Mq ( x) 


xH(x) 


(2.30) 


where H(x) is an analytic function of x. Then using Eg. (2.30) 
and defining a coefficient matrix by 


B(x) 


M 

_a 

H 


-V 


-s 


s 


(2.31) 
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the first of Eqs. (2.29) is rewritten as 


W = - BW . (2.32) 

X 

Next one looks for a solution of the system (2.32) in the 
Frobenius form 

CO 

W = x^ I C.x’^ , (2.33) 

k=0 


where X is an unknown constant and the are unknown 2x2 
constant matrices with Cq / 0. Since B(x) is an analytic 

function of x, it can be written as 


B(x) 


CO 


I 

k=0 


B, X 

k 


(2.34) 


where the constant matrices Bj^ can be found by expanding each 

of the elements of B(x) given in Eq. (2.31) in power series. 
Thus, for example, upon noting that V(0) = T(0), 


B 


0 


Mq (0) 

H(0) 


T(0) 

-T(0) 


-S(O)’' 

S(0) 


(2.35) 


Substituting expressions (2.33) and (2.34) into Eqs. (2.32) 
and equating corresponding powers of x yields 


(Bq - XI)Cq = 0 (2.36) 


and the recurrence relation 

k-1 

(B - (X+k)I)C = - I B C . (2.37) 

0 ^ j=0 ^ 


Since Cq 0, Eq. (2.36) implies that X must be an eigenvalue 
of the matrix Bq. By using the definitions of T and S from 
Eq. (2.26) and the fact that 
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1 - M^(x) 

liiti H(x) = lira [ ] = -2M' (0) 

x-^0 x->0 ^ 


the eigenvalues of are easily found to be 

= 0 


A2 = -1 - 


ifi[iX±JiI] 


( 0 ) 


1/2 


(2.38) 


Since and are distinct and do not differ by an integer, 
then corresponding to each A the eigenvector Cq can be deter- 
mined from Eq. (2.36) and the other follow from the recur- 
rence relation (2.37). In addition, standard theorems (Ref. 
16) insure that the resulting infinite series converge and 
are actually solutions of the system (2.23). Thus, two 
linearly independent solutions, and hence a general solution, 
can be found for the system. 

To determine the form of the general solution for x < 0 
it is assumed that M_ (x) is given by the extension (2.27). 

^ 1/2 

Then Mq( 0) = [(Y+l)a/2] , and two linearly independent 

solutions of the system are the analytic solution W^(x) 
corresponding to A^^ = 0, and, corresponding to 

1/2 

W“(x) = ^W^(x) . (2.39) 


Here the matrix W 2 (x) is analytic at x = 0 and W2(0) / 0. 

Similarly for x > 0 the extension given by Eq. (2.28) is 
used to obtain two linearly independent solutions: Namely, 

W^(x), which is analytic, and 


1/2 

TT+/ \ -(1-infi/a ' )A+/„\ 

W 2 (x) = X V ?2 (x) 


(2.40) 


Thus the general solution to system (2.32) can be written as 
the linear combination 
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(2.41) 


+ 

W = 


On 


On 


[ 


In 

— w;(x)e 
xl 2 


+inln X 

'1/2 


+ B£^W£ (X) 


where and are arbitrary complex constants, and the + 

and - signs are taken for x > 0 and x < 0, respectively. In 
order to satisfy general acoustic boundary conditions 

will not vanish, and thus the amplitude of the acoustic 

quantities will grow like x ^ as x 0. In addition, their 
phases are logarithmically infinite at x = 0 . One can observe 

that, even if A^^ = 0 jump discontinuities in the derivatives 

of r„ and at x = 0 will exist in general because of the 

On On 

jump discontinuity in (x) . 

In the following development, it will be required to have 

an explicit representation of two terms of the power series 

+ _ + 
for W“ and one term of the series for W^. Thus 



(2.42) 


where unified expressions for and have been introduced 

by use of the absolute value. The coefficient matrices in 
these expressions can be determined by solving Eqs. (2.36) 
and (2.37) for the x < 0 case and the analogous equations for 
X > 0. Equivalently both cases can be treated simultaneously 
by solving the matrix equations 
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9 


(Bq - XI) Cq = 0 
[Bq - (X+1)I]C* = -B^Cq 


(2.43) 


In Eqs. (2.43), X = -1 - (in^2 sgnx)/a^^^, B^ = B(0), and 

B^ = B' (0) where B(x) is given in Eq. (2.31). 

It is a matter of straightforward but lengthy algebra 
to compute the coefficients Bq and B^ and to solve the algebraic 

Eqs. (2.43) for the three required matrices. An outline of this 
computation is provided in Appendix A. The results of interest 
come out to be: 


•=02 = 




, (Y+Da/''^ 

- [ — ] sgn x 


(2.44) 


+ + 

, inS^ 

"b (3 -y) 

r n > 
2a 1 

1/2 -n 

^12 "^11 = 

1 + ^ 

_a 6 

[y+iJ 

sgn X 


,± ^ inS7(y + l) - 2 (Y - l)a^^^ sgn X 


‘02 


inQ (y + 1) - 4a^^^ sgn x 


in which a and b are the coefficients in the expansion (2.15) 
for the duct area ratio. 

Using the results of Eqs. (2.41), (2.22) and (2.24) the 

local form of the acoustic solution for any duct described by 
Eq. (2.15) can be written as 


r = A, 
n In 


= A 


X 11 


exp inS 




+ B^^ exp inQt + o (e) 


n In I X 


+ c^2 


( 4 . , sgn X 

r II7T 


In X 


(2.45) 


+ ^in^o 2 inJ^tt o(e) , 
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i i + 

where the values of ^02 given in Eqs. (2.44). 

It is the local expansion Eq. (2.45) which provides the 
key to the analysis of the acoustic behavior of the converging- 
diverging duct. For any fixed 6, no matter how small, the 
linearized acoustic theory fails to approximate the perturbation 
quantities near the throat of the duct as the flow Mach number 
there approaches unity. In general, both p and u become 
arbitrarily large in this circumstance, thereby violating the 
assumption made in deriving Eqs. (2.23) that they remain small. 
Numerical solutions of the lienarized acoustic equations 
illustrating the build-up of these arbitrarily large amplitudes 
were presented in Ref. 6. 

The expansion (2.19) will be termed the outer expansion 
in what follows. It is assumed to be asymptotically valid as 
an approximation to the solution to Eq. (2.4) as 6 -> 0 with 
X, t and e fixed. The result (2.45) shows that it is not 
uniformly valid as x ^ 0 and e 0. In the vicinity of x = 0 
nonlinear effects become important. The correct nonlinear 
theory which must be applied to describe the sound field in 
the region of near sonic flow near the throat will be presented 
in the next chapter. 
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IV. NONLINEAR THEORY 


A. Inner Expansion and Basic Equations 

The results of the previous section lead to the con- 
clusion that linearized acoustic theory fails to properly 
describe sound propagation in the vicinity of the throat 
for e 0. Both the singular nature of the outer expan- 
sion of the basic steady flow Eq. (2.14) and of the 
acoustic solution Eq. (2.45) indicate that the size of 
the region around x = 0 in which the outer expansions 
fail is of 0(e). Thus we introduce a "stretched" or 
inner variable 


X = I (3.1 

and examine the basic flow and the perturbation solution 
in the inner region, X = 0(1). 

Equations (2.2) can be used to derive the well-known 
equation for the flow Mach number in terms of the area 
ratio of the duct: 


^ - mg 1 ^ 

dx 1-M^ a dx ‘ ) 

The solution to this equation is given by Eq. (2.10), but 
it is convenient here to solve Eq. (3.2) over again in the 
inner region. An inner expansion for M of the form 


M(x;e) = M^(X;e) = m^ (X) + em^^ (X) + ... (3.3) 

is assumed. Transforming Eq. (3.2) to the independent 
variable X of Eq. (3.1), substituting the expansion (3.3) 
into Eqs. (2.5) and (3.2) and using Eq. (2.15) for a, then 
yields the sequence of equations on the m. . The first two 
of these are ^ 

dm^ dm 

■dir = “ ' '"i d5T ' • 

Thus, 

m^ = bQ , and m^ = ± (aX^ + b^)^'^^ , 


where b^ and b^ are constants which in general would be 
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determined by asymptotic matching of the expansion (3.3) 
to the outer expansion (2.18). However, in this case the 
matching process can be avoided because of the fact that 

the small parameter e is defined in terms of itself. 

Thus, M^(0;e) e 1- e, which implies that b^ = 1, b^^ = 1, 

and that the minus sign on m^ is appropriate. Hence, the 
inner expansion of the basic"^flow Mach number becomes 

M^(X;e) = 1 - e(aX^ + 1)^/^ + O(e^) . (3.4) 


Using Eq. (3.4), corresponding inner expansions of any of 
the basic flow quantities can be constructed. For example 
from Eq. (2.5) 


and 


G(X;e) = G^(X;e) = e(Y+l)m^(X) + O(e^) (3.5) 

M' (x;e) = ^ II = -(I±l)^/2 


where 


m^ (X) = (aX^ + 1)^/^ 


(3.6) 


Similarly 


U^(X;e) = c(X)M^(X;e) = c 


(X; e) 


[G^(X;e) 


. (3.7) 


The above results for the inner expansion of the basic 
steady flow will now be used to derive the non- 
linear inner equations for acoustic propagation. The 
Riemann Invariant form of the perturbation equation will be 
treated in detail while analogous results for density and 
velocity will be sketched in Appendix B, 

If the relations (3.4) and (3.5) and the change of 
variable (3.1) are substituted into Eqs. (2.7), the 
equations become 

T-1 

Pj + + ...)[ (1-em^ + . . . ) (1+u^) + (1+ P^) ^ ]PJ (3.8a) 

2in * 

" 7+r + •••] [(1+u"*-)^ - d+p"^)^"^] 

+ mjl^ (1 + em^ + . . . ) u^ [ (1-em^ + . . . ) (1+u^) + (1+p^) ^ ] = 0 , 
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and 


Y-1 


(|+^mi+ ...) [(l-emi+...) (l+u^)-(l+p'^) ^ ] Q J (3.8b) 


2m 


Y+1 


Y [1 + em^ + . . . ] [ d+u"-) 2 - (1+p^) 


Y-1 

+ (1 + ■*■ • • • ) [ (l-em^^ + . . . ) (1+u^) - (1+p^) ^ ] = 0 


respectively. In Eqs. (3.8) the notation u^, p , P and 
has been used to denote the expressions for the perturbation 
quantities as functions of the inner variables X and t. Now 

it is assumed that u^ and p^ can be expanded in powers of e 
as 

u^(x,t;e) = ey^(x,t) + e^y 2 (x,t) 
p^(x,t;e) = er^(X,t) + e^r^CX/t) 


(3.9) 

t . . . . 


These expansions of u and p yield the corresponding expres- 
sions for P and Q from Eq. (2.6): 


where 


"'(X, 

t; e) 

= eP^ 

(x,t) 

+ 

(X,t) + 

/ 







(3.10) 

'^(X, 

t ; e) 

= eQjL 

(x,t) 

+ 

(X,t) + 

• 


i 

i 


^ i 

i 


hd 

H 

II 

^1 


f 


- ^1 ' 


^2 = 

i 

^2 

- 4 - 

Y-3 

4 

(r^)2 - 

i 

m^^ y f 

(3.11) 


i 

i 

y-3 

/ i\ 2 

i 


Q2 = 

^2 

^2 

4 

(rp) - 

m^y^ 



The expansion for P from Eq. (3.10) and the series of 
Eq. (3.9) are substituted into Eq. (3.8a). Equating the 
coefficients of powers of e separately equal to zero then 
yields 
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0 


(3.12a) 


0(e): 


0(1): = 



Y-1 

y+1 


m|Pi 


(3.12b) 


Similarly, if the expansion for from Eq. (3.10) and the 
series of Eq. (3.9) are substituted into Eq. (3.8b), the 
lowest order terms yield 


+ ( 


3-y 


Y+1 

4 


Qi- 


m^) 


X 


2mj^ 

Y+1 


■3-y 




Y+1 

2 


Qll = 


0. (3.13) 


The inner equation of motion on follows in a similar manner. 

However, because it is complicated and will not be needed in 
the following, it is omitted here. 

Equation (3.12a) and Eq. (3.13) consitute the leading 
order nonlinear theory governing sound propagation through a 
near-sonic-throat of the form of Eq. (2.15) with a ^ 0. At 
this point, considerable progress in treating these equations 
can be made by taking advantage of certain facts which are 
apparent. First, both the first and second order equations 
on P are linear. Second, the combination y + r from Eq. (2.45) 

in non-infinite as x 0, because the terms in |x| ^ cancel, 
whereas y - r in the outer (linear acoustic) solution is 
infinite as x ->■ 0 . These considerations suggest that in 
reality P is of smaller order than Q in the inner region, 
and that the downstream wave passes through the throat essen- 
tially unaffected by nonlinear interactions with the near- 
sonic flow. This suggestion is supported by physical reason- 
ing as well. The Q wave, which propagates against the flow, 
is intensified relative to the downstream wave P by a 
transonic trapping near the throat. Thus, it will be assumed 
here, and verified later by asymptotic matching, that P^ s 0 . 

Then Eqs. (3.12) and (3.13) imply 


Y+1 






and that P 2 is independent of x. 
If the definition 


(3.14) 


TT(X,t) 


_ Y+1 
4 


Qn 


m. 


(3.15) 
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2 1/2 

where = (aX +1) f is introduced into Eq. (3.13) it 
becomes 


Tr^ + = aX . (3.16) 

The solution of this equation will be discussed in the next 
section. 


B. Solution of the Nonlinear Inner Equation 

The nonlinear partial differential equation 


IT + 1T1T = aX , -oo < X < oo (3.17) 

t X 


which describes the upstream propagating sound in the throat 
region will now be solved by the method of characteristics 
Ref . 17„ It follows directly by this method that the 

characteristic curves in the t, X plane are solutions to 
the differential equation 


dX IT 


and that along these curves. 


dir 

dX 


' ^ ' 

r dt ■ 


l dx ^ 


aX 

TT 


(3.18) 


Integrating this ordinary differential equation along the 
characteristics yields 


TT = ± (aX + J) 


1/2 


where J is a constant along a characteristic. It can be 
shown using matching considerations similar to those given 
in Section III-C below that the solution to the above 
equation corresponding to the plus sign must be discarded. 
Thus for the remainder of this report the solution to be 
used is 

IT = - (aX^ + J) . (3.19) 
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If positive solutions for ir are assumed then it can be 
shown that the total axial velocity u will be supersonic 
throughout the inner region. Hence, unless shock waves are 
introduced in the inner region, matching with the subsonic 
outer solution is not possible. [It is expected that positive 
solutions will be necessary to treat cases with shocks present.] 

It is convenient in solving the specific problems con- 
sidered in Section III-D to introduce the parameter t, where 
t = T along the line X = 0 and to call the value of it at X= 0, 

4) (t) . Then using 


TT(0,t) = - = (f)(t) = 4>(t) 

the solution tt can be written in the parametric form 

u(X,T(X,t)) =- [ax^ + 4)^ (T) ] . 


(3.20) 


Using Eq. 
separable 
parameter 


(3.20) in Eq. (3.18) and integrating the resulting 
ordinary differential equation yields the one 
family of charasteristic curves 


t = T - 


1/2 


In 


1 / 2 ^ ^ r ^2 ^ ^ 2 , , , 1/2 

a ' X + [aX +4) (t ) ] ^ 

I 4) (t) I 


(3.21) 


If the function (p is known then Eqs. (3.20) and (3.21) 
provide a solution to Eq. (3.17) for -t» < X < «>. This will 
be the case for an acoustic source of 0(e) located at the 
throat of the duct; hence, the inner solution is completely 
determined to leading order and plays the role of a boundary 
excitation on the acoustic field in the outer region. On 
the other hand, if the source of sound is located in the 
outer region, then 4>(t) cannot be determined without 
simultaneously solving the outer equations and asymptotically 
matching their solutions with the inner solution tt. Regard- 
less of the location of the source (p (t) must be a periodic 
function, with period 2 tt/^^ equal to that of the source, since 
the overall acoustic field and hence tt must be periodic 
functions. In the remainder of this section and in the next 
section on asymptotic matching 4> shall be assumed to be a 
periodic function. When 4> is periodic then tt(X,t) is periodic 
in T with the same period. In addition Eq. (3.21) implies 
that t - T is a periodic function of t with this period. Thus 
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TT is periodic in t with the period of cj>. 

Before carrying out the matching several observations 
are in order concerning the solution tt and the characteristic 
curve. In particular a condition will be derived which 
insures the existence or nonexistence of acoustic shocks. 

If IT is allowed to vanish, then a careful analysis would 
show that 8 it/ 3X will be infinite, and hence shocks in the 
solution TT are to be expected. In addition it follows from 
Eq. (3.20) that 


9tt 

3X 


3it 

3X 


+ 


where 


3ir 

3t 


X 


4)' - 


3t 

3X 


d({) 

dx 


aX 4>c}) ' 3 t 

TT ■ TT 3X 


(3.22) 


Rewriting Eq. (3.21) as 


F(X,t,7T) = t - T - 


1/2 


In 


1/2^ ^ / v2 ^ a2^ 1/2 

a X + (aX + (j) ) 


= 0 


and using the implicit function theorem yields 

3t _ “^X _ (|) 

3X TT(j) - X(f) ' 

where 


F 


X 


1 

TT 


and 


F = - 1 + ^4)' 

T TTCt) 


(3.23) 


(3.24) 


It is clear from Eq. (3.22) that if 3x/3X is finite then 
3 tt/ 3X will remain finite (and hence no shocks will appear 
in the inner solution). It follows from Eq. (3.23) that 
3x/3X remains finite if 


Sq = TT(j) - X(f) ' 


2 2 

- (aX +<j) ) 


4> - X(t>' ^ 0 


(3.25) 


Since at X = 0 Sq = - 1 4> | 4> > 0 / Sq> 0 for all X _< 0 , t > 0 . 
Using Eg. (3.25) the condition S^> 0 is equivalent to 


- ^.2 ^ ,2,1/2 

(aX + (j) ) ^ 


Xcp ' 

4 > 
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This condition is always satisfied if <j)' >0 since < 0. 

If however, (j) ' < 0 for any t which is the case for c() (t) 
periodic then > 0 if and only if 

acj)^ - ((f)')^ ^ 0 . (3.26) 


If condition (3.26) is not satisfied, then it follows 
from Eq. (3.24) the vanishes at some X,t. Hence, the 

characteristic curves intersect. Therefore, violation of 
condition (3.26) implies that shocks must appear in the 
inner perturbation quantities. The physical interpretation 
of Eq. (3.26) will be discussed in Section III-D. 

C. Asymptotic Matching 

The solution in the inner region has been expressed in 
a form suitable for matching in the previous section. This 
solution will now be matched to the outer solution in order 
to relate the overall perturbation level 6 away from the 
throat region to the perturbation level e in the throat or 
inner region and to relate the acoustic field in the outer 
region to the acoustic field near the throat. The matching 
process will thus involve a distinguished limit and is 
intricate. Since the basic steady flow quantities have 
been matched in Section IIIB, it is only necessary to match 
the perturbation velocity and density u and p or equivalently, 
the perturbation Riemann invariants P and Q. 

Here the matching will be carried out for P and Q. The 
matching of u and p is assured by the matching of P and Q. 

For details of the former see Ref. 4. In order to completely 
determine the leading order inner and outer solutions matching 
of the leading order P and Q waves will be required. 

The one-term inner expansion of Q, denoted here by Q , 
is given by Eqs. (3.10), (3.15) and (3.20) in the form 


= sQ^(X,t) = ^ [(aX^ + 1)^/^- (aX^ + 4 >^(t))^/^] . (3.27) 

To apply the matching principle of Van Dyke, Ref. 18 , 
to relate the inner to the outer solutions a one-term outer 
expansion of the expression (3.27) is required. This is 
obtained by using Eq. (3.1) to express the inner quantity 
in outer variables 
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7+T X 


t 


4 

"a (y+1) " 

1/2 r 

1 V 1 

7 2s2 ] 

1/2 

f „ 2,2 

2e (j) 

1/2- 

Y+1 

2 

1 X 1 

^ ^ 2 
a(Y+l)x / 


^ 2 
a(Y+l)x 

- 


Expanding this expression for fixed x and small e using the 
series 

(1 + y)^/^ = 1 + |y + ... 
and retaining the leading term in e yields 



, 2 3/2 



[1 - (T) ] 


(3.28) 


Here the notation is used to indicate the one-term outer 

ol 

expansion of the one-term inner expansion . 

The parameter t in Eq. (3.28) is determined from the 
characteristics of Eq. (3.21), also by expanding in e with x 
fixed and retaining only the leading term: 


where 


g 


t = T + ln|(|)(T) I + g(x;e) 


(x,-e) ln[ 2 ( 15 i a) 1/2 M, 


(3.29) 


To utilize the results (3.28) and (3.29) requires dis- 
playing the explicit dependence of on x and t. Because 

cj) (t) is a periodic function of t with period 2ir/n it can be 
expanded in a Fourier series as 


♦ (T) = jC- + 


n=l'- 


+ + 

C (x) cos nOt + D (x) sin nfit 
n ' n ' ' 


(3.30) 
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in which the + and - signs refer to x > 0 and x< 0, respectively. 
The evaluation of the coefficients in Eq. (3.30) is carried out 
in detail in Appendix C. The results are 


n^fl^ 1 1 * I cos nf^ (g±e*) 


C- = ^ (1+ii^) 
n 7T 4a 


2 2 - 1/2 

D" = - ( 1 + [i*! cos nn(g±ej) 
n IT 4a ' n‘ n 


In Eq. (3.31) the are defined by 


(3.31) 


+ 

I" = 
n 


2i\/n 


(j) I ^ exp (±inJ^x) dx 


(3.32) 


and the phase angles by 


+ 

H“ 

1 “1 nfl n 1 

0 o [tan =r-PT t ~T J 

n nJ^ T?~ 

2a 


(3.33) 


n 


+ + 

in which E" and H“ are the real and imaginary parts of I 
n n 

The exponent v in Eq. (3.32) is 

infl 


V = 2 - 


1/2 


Therefore, combining Eqs. (3.28), (3.29), (3.30), and (3.31) 

the explicit representation of is 


Q 


il 

ol 


, 2 , 3/2 
>+l^ 




[1 — ^ 1 1~ 
2tt I-^O 


00 

r 

' ^ nil 


2„2 - 1/2 
(1 + ^^^) ’ 


1^1 cosnfi(t-g + 9 )]. (3.34) 


The outer expansion of the complex solution Q for a 
periodic wave of period 2ir/n can be constructed as a Fourier 
series from Eqs. (2.21), (2.22) and (2.24) in the form 
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CO 

Q° = (S ^ exp (inf^t) +0(e6) . . (3.35) 

n=l 


The one-term inner expansion of the leading term of this 

outer expansion, denoted here by is formed by expanding 

each harmonic component for small e, keeping X fixed, and 
retaining only the leading term. This is equivalent to 
expanding each harmonic for small x, which is precisely the 

result given for y^^ and r^^ by Eq. (2.45). If Eq. (2.17) 

is used for (x) , then to leading order for small x Eq. (2.45) 
yields 

-2A* 

Von-'^On - • 

' ' a 

The physical value of Q° is found by taking the real part of 
Eq. (3.35). Thus, if n ^ 1, is written as A“ exp (inJ^'i'~) 

with A“ real, the one-term inner expansion of the physical Q 
follows as 


^il ^^10^ I A%osnfi(l+ ln|x] +¥^)] .(3. 

' ' n=l a 
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Now, by the matching principle. Ref. 18 = Q°^. 

Equating the results of Eqs . (3.34) and (3.36) gives 

X 2 

0 = e 


A± =-l (2)3/2 

10 2 V+1 


A“ - — 

n 2 tt 'y+1 


(r?^) 


1/2 2tt 


. 1/2 


2^2 -1/2 

(1 + ^1/^) 

4a 


n' 


(3.37) 


’^n }~T 72 


, r 2 ,Y+1 ,1/2- 

^ e ^^2 ^ ~ ®n| ^ 


Eqs. (3.37) yield the relationships between the complex constants 
+ 

Aj^^ of Eq. (2.45) in the linearized solution on either side of 
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the throat region. In addition, the result shows that the up- 
stream wave Q, which is 0(6) in the outer region, is strength- 

1/2 

ened to 0(e) = 0(6 ) in the inner region as a result of non- 

linear interactions with the near-sonic flow. 

The results of Eqs. (3.37) indicate types of nonlinear 
effects predicted by the theory. They include in addition to 
the intensification of the Q-wave a generation of superharmonics 
and an acoustic streaming effect. 

It remains to determine the relation between the complex 
+ 

constants of Eq. (2.45), which is done by asymptotically 

matching the downstream wave P. In the inner region the fact 
that s 0, and Eq. (3.12b) imply that P 2 = ^2 thus 

using Eq. (3.10) implies 


P""^ = e^P2(t) 

which is equal to P^^ because P 2 is independent of X. On the 
other hand, the complex P in the outer region is 


P°1 = S I (MQPon ^On^ (inQ) 

n=l 


(3.38) 


from Eq. (2.21). Use of Eq. (2.17) for Mq(x) and expansion of 

the components of the series (3.38) for small x using Eq. (2.45) 
yields for the leading term 


^0^0n'*'^0n = ^In ^^11 ^12 ^^2 ^ ^ sgn x] 

X exp infl (t + ln|x|) + 




Then Eq. (2.44) shows that the term multiplying A, vanishes 
in general so that 



00 

6 I B" (1 + d* ) exp (inJ^t) 
n=o 


(3.39) 


in which d ^2 is given in Eq. (2.44). This result verifies the 
earlier assumption that P^(X,t) = 0, as the leading order outer 
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term for P is not singular as x ->■ 0. Finally, since P^^ is 

independent of X, the matching can be accomplished by simply 
equating the complex expressions (3.39) on either side of 
X = 0 . Thus 






and solving for the ratio of to using Eq. (2.44) for 

dp yields 


^In _ [in^^(Y+l) + 4a^^^] (in^^ - a^^^) 
bI!' [inf^(Y-l) - 4a^^^] (inf^ + a^'^^) 


(3.40) 


Equations (3.37) and (3.40) represent the primary results 
of the inner solution process. If the function 4> (t) is known, 
then they provide two complex equations for the four complex 
+ + 

constant and B^^ for each n. The remaining two equations 

to determine these constants are the acoustic boundary and 
source conditions which will be stated for specific cases in 
the following section. In general, of course, (j) (t) is not 

+ 

known, and must be determined along with the constants and 

B^^ taking account of the complicated coupling between the 

inner and the outer solutions. The general case can ^e reduced 
to a coupled system of transcendental equations for and 'i'^. 

It will not be discussed further in this report. 

D. Solution of Specific Boundary Value Problems 

In order to complete the specification of a boundary 
value problem two boundary conditions must be prescribed for 
the perturbation quantities. These could both be applied in 
the outer region, on either side of the throat. A typical case 
would be an anechoic termination upstream of the throat and a 
single harmonic sinusoidal velocity source downstream of the 
throat. Another possibility is to have one boundary condition 
applied at the throat of the duct and the other either upstream 
or downstream of the throat. Several examples of the latter 
case will be discussed in this section. 

As was discussed earlier, the function <)) (t) , upon which 
the asymptotic matching is based is known a priori only if the 
acoustic source is located at the throat of the duct. In order 
to illustrate the theory in the present report, it will be 
assumed that cp is known. This allows solution of two specific 


36 



problems; transmission of sound upstream or downstream from 
a source at the throat and transmission through the throat 
from a reconstructed source in the outer region which yields 
the specified <p at the throat. 

In order to solve these problems it will be assumed here 
that the function cj) (t) is given by 


4>(t) = (1 -Nq)(- 1 + N sin nt) (3.41) 


where Nq and N are specified constants. It will now be 
verifiea that this 4> corresponds to a perturbation velocity 
at the throat having both a steady and a single unsteady 
sinusoidal component at frequency fl. To show this consider 
a source condition on the total perturbed velocity field of 
the form 


U(x;e)u(x,t;e) 


x=0 


(X; e)u^ (X,t; e) 


X=0 


/V. y\ 


Uq + U sin t^t 


(3.42) 


/N /N 

where U and Uq are constants. Here the source strength is 

taken to be of order e to be consistent with the condition 
that the acoustic perturbation quantities are of order e in 
the inner region. Using Eqs. (3.7) and (3.9) this condition 
becomes to first order in e 


y^(0,t) = ^ (1^)^/^ (Uq + U sin f^t) (3.43) 

s 

since ~ + r^ = 0 and ~ = 2y^, Eq. (3.43) 

implies 

Q^(0,t) = (UQ+Usinfit) 

and from Eq. (3.15) 

iT(0,t) = Q^(0,t) - m^(0) 

= ^ (^^) ^'^^ (Uq + U sin J^t)-1 . (3.44) 
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I II 


mill II 


However iT(0,t) = (t) when (j> (t) < 0 and then equating Eqs. 

(3.44) and (3.41) gives 


N„ = f— 

0 c 




(3.45) 


Cs(l-N„) 


(3.46) 


Thus prescribing a velocity source of the form Eq. (3.42) at 
the throat of the duct completely determines c}) (t) in the form 
of Eq. (3.41). It should be noted that if = 0, the applied 

source has no steady part. 

Since (() (t) is known, use of Eq. (3.41) in Eq. (3.32) gives 
So that Eq. (3.37) yield for A“ 


+ 12 1 2 

cl 


(3.47) 


Now the coefficient A^q represents a steady component in the 

acoustic perturbation quantities in the outer region, which is 
a form of acoustic streaming. If N = 0 then there is no such 

^ ± 

steady component at the throat. In this case A. „ 7 ^ 0. However, 
by choosing 


,2 - 1/2 


Nq = 1- (1 + ^) 


(3.48) 


the streaming term A^q can be suppressed in the outer region. 

One case to be considered here will require that condition 
(3.48) holds with / 0. Thus no streaming will occur in the 

outer region. The other case to be considered in detail will 
have a nonzero streaming term. 
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If Nq = 0 then A^q ^ 0 for <j) (t) given by Eq. (3.41) and 

a streaming term will occur in the outer solution. Here an 
acoustic source at the throat of the duct will cause an acoustic 
stream effect in the outer region. 

Since solutions without shocks are being considered in the 
report the function (j) (t) must satisfy condition (3.26) namely 

a((|))^ - (<|)*)^ > 0 


If (j) is given by Eg. (3.41), then, this condition becomes 

a . cos fit 2 

„2„2 — '■1 - N sin flt"^ 
iz N 

Maximizing the right side of this inequality with respect to t 
yields 

^ , or < P (3.49) 

N 1-N 1 + ^ 

a 

Thus in order to have no shocks in the inner region the 
conditions <j) < 0 (since tt < 0) and Eq. (3.49) must be satisfied. 
The condition (j) < 0 implies that if = 0 then N '< 1 while if 

Nq 0 then 0 < Nq < 1 and 0 < N < 1. It is interesting to note 

in the case Nq = 0 that if the parameter e is reintroduced then 

the source strength is in effect eN. If this product is denoted 
by Sq then Eq. (3.49) can be rewritten as 



a 


In this form the shocking condition relates the source strength 
and frequency, the curvature of the throat and the deviation 
of the throat Mach number from unity. As is to be expected 
shocks will occur at lower source strengths for higher Mach 
number flows. Similarly a high frequency source or a flat 
throat region (small a) will lead to earlier shocking. 

To completely prescribe the boundary value problem and 
+ + 

determine A^^ and two more conditions are required. These 

are the acoustic boundary and source conditions. If the 
boundary conditions are applied well away from x = 0 then it is 
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necessary to solve the outer equations (2.25) before completing 
the solution. This can be done numerically. To illustrate the 
theory here, however, it is more convenient to choose a special 
dust shape due to Crocco and Tsien for which an analytical 
solution to the system (2.25) is known. This solution is 
discussed in Appendix D. 

For the case of upstream propagation from the throat it 

will be assumed that the plane x = - x of the Crocco-Tsien duct 

e 

is anechoic so that the acoustic boundary condition there is 


P°^(-x ,t) = 0 . (3.51) 

e 

Use of Eqs. (3.38) in Eq. (3.51) as discussed in Appendix D 
then gives 


E 


— = (2K)^"^ = r (-X ) 

n e 


2n 


B 


In 


(3.52) 


in which T^(x) is a ratio of combinations of hypergeometric 

functions whose specific form is given in Appendix D. 

Eq. (3.52) is the third relation used in the determination 
+ + 

of A^^ and in the case of anechoic conditions upstream 

of the throat. ^ ^ 

The final equation to determine A“^ and is the source 

condition. For the source at the throat, prescribing the 
acoustic perturbation velocity results in Eq. (3.41) for (}) (t) , 
which completely solves the problem of proagation out through 

the exit upstream. A^^ is given by Eq. (3.37) directly, 

is determined using Eq. (3.47), and the sound field is 
constructed using Eqs. (d 7) and (D8) . 

A similar procedure would solve the problem of downstream 
propagation from a throat source if the exit condition were 
applied at a positive x. In this case Eq. (Dll) is used ^ 
instead of Eq. (3.47) and the source condition determines A^^ 

instead of 

The procedure just described was implemented numerically 
and results were obtained for three different types of boundary 
value problems. The details of the numerics are discussed in 
the next section. 
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For the special case of transmission from a source down- 
stream of the throat which can be treated using the above 
approach, (t) is still assumed to be given by Eq. (3.41). In 

this case Eq. (3.37) yields Eq. (3.52) gives and 

Eq. (3.40) determines The acoustic field on either side 

of the throat is then given by Eqs. (D7) and (D8) . This result 
can be viewed as the sound field transmitted upstream from a 
periodic velocity source located at some x > 0, where the 
strengths of each of the harmonic source components are 
constructed using Eq. (D8) at the source location. This pro- 
cedure is special in that the source strength and harmonic 
content is not to be specified but rather is part of the 
solution. An iterative procedure could be set up to control 
the harmonic content of the source by changing the form of cj) (x) 
(by increasing its harmonic content) . This case will not be 
discussed further in the present report. 

E. Numerical Results and Discussion 

The results to be presented here were obtained by 
numerically implementing the procedure described in the last 
section. That is: 

(1) A choice of the parameters N^, and in <j) (t) 

(Eq. (3.41)) and the parameter a was made. Then the integrals 
were evaluated numerically by Simpson's rule. This deter- 
mines, using Eq. (3.37) , the constants A^q, and upon choosing 
a value of e the phase angles 

n 

(2) Either Eq. (3.47) or Eq. (Dll) was then employed to 

determine either or B~^ respectively. The former being 

used when the anechoic condition is applied upstream of the 
throat and the latter for a downstream anechoic boundary. The 
hypergeometric functions in these expressions were evaluated 
by use of the standard power series representation for these 
functions . 

First although the results of the asymptotic matching are 
sufficient to transfer the acoustic field through the duct 
without further consideration of the inner solution, it is of 
interest to display this solution to illustrate the approach 
to shock formation discussed earlier. Figures 5-8 show 

the inner acoustic velocity y^(X,t) evaluated by use of Eqs. 

(B7) and (3.21) at one fixed time for |x I £ 10 and for the 
indicated values of the parameters N, N^, and a. In Figure 

5 the value of = 0 and hence the source at x = 0 has no 


41 


steady part. The steepening to a shock is clearly indicated 
and condition (3.49) yields the value N = .7746 for shock for- 
mation. Figure 6 yield similar results for the indicated 
values of N, with related to N by Eq. (3.48). Similarly in 

Figures 7 and 8 the variation of y^(x,t) with the source 
frequency Q and the geometric parameter a respectively is 
shown. In both cases the shock formation is evident. 

In Figure 9 the acoustic velocity is shown for six 
equally spaced times over a period 2-n/U. The upstream propa- 
gating nature of the inner acoustic velocity is clearly evident. 


The remaining discussion will require an expression for 
the time averaged acoustic power transmitted out through the 
exit of the duct which in each case is taken at either x = +.06 
or X = -.06. The acoustic power, calculated according to the 
expression given in Ref. 19 can be written in terms of the 
dimensionless variables of this report as 


W 

3 

R c 
s s 


a (x)Mq 
2 


fl - 3t] 

r- 


[2 (Y-1) J 

2 

[^^00^ ‘‘‘^O^^OO^ ^^"^^0^ ^^00 ^ ^^00 


+ 



+ m: 



+ 


(1+Mq) 




cosnf2(6 -0 

r n y nj 


in which 0 and 0 are the phase angles of the complex r„ and 
y^^, respectively, and R^ is the stagnation density of the 
steady flow. Here r^^ and y^^ correspond to steady perturba- 
tions or acoustic streaming type terms. 

Figures 10 - 13 give results for a source located at the 
throat of the duct, with = 0 and an anechoic upstream 

boundary. In Figure (10) results for two cases with the source 
strength fixed at = eN = .03 and .06 are presented. The 

effect of increasing the throat Mach number (decreasing e) on 
the upstream transmitted acoustic power is determined. In both 
cases there is marked attenuation with increasing Mach number. 

It is of interest to note that shocks will occur for e = .08484 
when Sq = .06 and for e = .04242 when = .03 (Eq. (3.50)). In 

Figures 11 and 12 the upstream transmitted power is pre- 
sented for varying and varying a respectively, with and e 

fixed in both cases. Results indicate that the transmitted 
acoustic power increases with increasing source frequency and 
decreases with increasing a. 
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In order to study the harmonic content of the solution in 

the outer region upstream of the throat the ratio W /W was 

n 

calculated where is the acoustic power carried by the n-th 

harmonic. In Figure 13 this ratio is plotted against the 
deviation of the throat Mach number from unity (e) for the 
fixed source strength = eN = .03. Curves are given for the 

first three harmonics. It is seen that as e decreases the non- 
linear interaction causes an increasing fraction of the total 
power to be distributed among the higher harmonics. For 
e = .05 approximately 86% of the power is carried out in the 
fundamental harmonic. This nonlinear effect occurs even though 
no shock has formed in the inner region. It is quite 
important to emphasize that if linearized theory had been 
applied in the inner region all of the acoustic power would have 
remained in the fundamental harmonic. This result is of obvious 
signficance in the design of experiments to measure sound radia- 
tion out of near-sonic inlets. 

The remaining results discussed in this report correspond 
to cases where N_ 0 and the parameters N and N_ are related by 
Eg. (3.48). Hence no streaming will occur in the outer region, 
in these cases. In Figures 14 - 17 results are presented for 
an upstream anechoic boundary. The results of Figure 14 are 
similar to those of Figure 13 in that about 15% of the energy 
is carried by higher harmonics for values of N below the value 
where shocks will occur. In Figure 15 , e = .1, = .01 and 

a = 1.0. Hence the throat Mach number and the source strength 
are fixed and the variation of W^/W with frequency is given for 

the first three harmonics. The graphs show that the nonlinear 
effect due to increasing source frequency leads to multiple 
harmonics and a spreading of part of the power to higher 
harmonic s . 

A similar result is shown for variable throat curvature 
parameter a in Figure (16) . In this case the flatter the 
throat the more pronounced the nonlinear effect. This is as 
should be expected since for a fixed throat Mach number 
(e = . 1 in this case) the duct with the flatter throat region 
will have a greater portion of high subsonic flow around the 
throat and hence a more pronounced nonlinear effect. Figure 17 
gives the dimensionless acoustic power W/R c^ radiated upstream 
as a function of for e fixed. As is to^be expected the 

acoustic power transmitted increases as the value of is 
increased . 
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In Figures 18 - 20 results are given for the case where 
the anechoic boundary was placed downstream of the throat and 
thus the power flow will be downstream. There is however an 
upstream propagating Q-wave in the duct because of reflections 
of the downstream propagating P-wave. These reflections are 
due to flow gradients (not the boundary) and occur in the 
outer region (downstream from the throat) . The presence of 
the Q-wave can be expected to lead to nonlinear results. 

This is in fact the case as is indicated in Figures 18 - 

20 where results similar to Figures 14 - 16 respectively 

are given. It is interesting to note that the total power is 

larger in the downstream case than in the upstream case for 

the same source at the throat. This is to be expected since 

in the former case the power flow is in the fluid flow 

direction. The spread of power to the higher harmonics is also 

more pronounced in this case. It is interesting to note in 

Figure 19 that as Q. increases the second and then the third 

harmonics have maximum values for W /W. This indicates as do 

n 

other numerical results that as shocking conditions are 
approached a great many harmonics are carrying the acoustic 
power . 
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V. CONCLUSIONS 


I. The study has shown, using a quasi-one dimensional 
model, that linear acoustic theory is not valid in general for 
sound propagation in near sonic duct flow. Linear acoustic 
quantities become infinite in such flows at the throat of a 
converging diverging duct. 

II. A new nonlinear theory to describe sound propagation 
in such situations is derived by the method of matched asymp- 
totic expansions. The nonlinear theory must be used in the 
vicinity of the throat while linear theory still holds away 
from the throat. Asymptotic matching is used to connect the 
nonlinear and linear regions and to completely determine the 
acoustic field in the duct. 

III. Analytical solutions are obtained to the nonlinear 
equations, thereby illustrating the physical propagation 
process. These are: 

(1) Intensification of upstream propagating sound due to 
a transonic trapping. 

(2) To leading order, downstream propagation obeys linear 
theory. 

(3) An infinite number of superharmonics are generated 
by a single frequency acoustic source. 

(4) Acoustic streaming effects are present in general. 

(5) Shocks will occur in the acoustic quantities when 
condition (3.26) or (3.50) is not satisfied. This condition 
relates source amplitude and frequency, steady flow Mach 
number and the curvature of the throat of the duct. 

(6) The acoustic power radiated out of the duct depends 
on the and thus decreases rapidly as the throat Mach number 
approaches unity. 

IV. Numerical studies are carried out for several simple 
boundary value problems to illustrate the predictions of the 
theory. These results include: 

(1) Appreciable percentages (over 15% in many cases) of 
the acoustic power transmitted upstream out of the duct will be 
carried by superharmonics of the source frequency. This can 
occur before acoustic shocks appear. 

(2) Increasing the source strength or frequency or de- 
creasing the geometric parameter a increases the nonlinear effect. 

(3) Shocks will occur at typical acoustic source fre- 
quencies and amplitudes and thus their effect must be studied 
further. 

V. Numerical and analytical results indicate that super- 
harmonics and shocks will be observed in near sonic duct flows. 
Experimental work to detect such effects in actual flows would 
be of much interest and benefit in understanding the acoustic 
processes at work. 
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Appendix A. LINEAR ACOUSTIC THEORY 


In this appendix a sketch of the derivation of the 
constants in the local expansion of linear acoustic theory 

Eqs . (2.44) will be given. To determine d“„ the last of 

Eqs. (2.43) 


“ 1 “ 

o 

o 

T(0) 

-s( 0 ) ' 


" 1 " 

1 

CM 

+1 o 

1 

H(0) 

-T(0) 

S(0) 


_ ^^02 _ 


0 (Al) 


must be solved. Here T and S were defined in Eqs. (2.26). 
This is just two linear homogeneous equations in two unknowns 
and thus eigenvectors of the matrix corresponding to the 

zero eigenvalue X = 0, are being found. Thus 


± _ T(0) 

02 S(0) 


1/2 ^^0 
^"''=0 -gT 


1/2 (Y 

inflGo + g 


X = 0 


(A2) 


Using the results 


G„(0) 


-I(Y+1) sgnx 

(Y+1) 

2 


(A3) 


Eq. (A2) reduces to the desired result (Eq. (2.44)). 

+ 

The constants are found by solving the first of 

Eqs . (2.43): 
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[B 


0 


- XI] 



0 


where 
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APPENDIX A 


X = -1 


(infl sgn x) 


(A4) 


and B_ is as given in Eq. (Al) . The solution is easily found 
to be Cq 2 = “1 tiy using the definitions of T(0), S(0) and Eqs . 
(A3) . 

To find c" and c“ the second of Eqs. (2.43) 



^C7, “1 


” 1 “ 


11 



[Bq - (X+DI] 

+ 1 
c, ^ i 

= -®1 

-1 


L 12^ 




(A5) 


must be solved. Here the eigenvalue A is given by Eg. (A4) , 
Bq is given in Eq. (Al) and hence the matrix B^ = B' (0) must 
be found. The computation of B^ is long but straightforward. 
From the definition of B(x) (Eq. (2.31)) it follows that 


(%] 

1 

~T 

-s' 

M 

+ ^ 

t' 

-t' 

H J 


j-V 

S _ 

^ IT 

_ -V 

s' 


or 


B' (0) = 


F'(0)T(0) + F(0)T'(0) -(S(O)F'(O) + S'(0)F(0)) 


-(F' (O)V(O) + F(0)V (0) ) S(0)F' (0) + S '(O)F(O) J 


(A6) 


where F(x) = Mq(x)/H(x). Using Eq. (A6) in Eq. (A5) and solving 
the latter by inverting the coefficient matrix yields 


+ 



c_ . 



11 

= r 

Aj^ + aA2 

+ 


* _ * 



B ~ aB^ 

_ 12 _ 


_ 2 


(A7) 


where 
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APPENDIX A 


a = ( 1 +X 2 )H( 0 ) , 


r = aH(0) [a - T(0) - S(0) ] 


A, = 


S(0)F(0) (V (O)-T' (0) ) , = T(0)F(0) (V (O)-T' (0) ) 


A 2 =F'(0)[T(0)+S(0)] +F(0)[T'(0)+S'(0)] 
B 2 = F'(0)[V(0)+S(0)] + F(0)[V'(0)+S'(0)] 


To arrive at the final result (Eqs. (2.44)) the quantities F, 

T, S, V and their derivatives must be evaluated at x = 0. This 
is tedious and the details will not be presented. It should be 
noted that it is easier to do the actual computations by 

+ + 

rewriting Eq. (A7) to obtain expressions for + c and 

± ± 

~ ^11 hence this is the form of the results presented 

in Eqs. (2.44) . 
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Appendix B. INNER EQUATIONS FOR VELOCITY AND DENSITY 

In order to derive the inner equations for the acoustic 
density and velocity Eqs. (3.11) could be solved for and r^^ 

in terms of and Q^ and the result substituted in Eqs. (3.12) 

and (3.13), or equivalently the expansions (3.9) and the steady 
flow expansions (3.4), (3.5) could be substituted into Eqs. 

(2.4) and coefficients of each power of e are set equal to 
zero. The latter procedure yields for the first of Eqs. (2.4) 



0(e): 
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0(e 

V i . i„i 

^ ' ^It ^1 IX 


for the 

second of Eqs. 

(2.4) 


0(e) : 

^^x 
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+ ^2X 


= 0 


X 

IX 


IX 

2aX 


+ [(y-2)r^+m^]r^^ - [2y^ 


(B2) 


(Y-1) r^] 


= 0 


The continuity and momentum equations yield redundant 
equations for the leading order density and velocity. Thus 
2 

the 0(e ) equations will be needed to derive two independent 

* * 2 
equations for r^ and y^. The redundancy is repeated at 0(e ), 

and thus subtraction of the second of Eqs. (Bl) from the 
second of Eqs. (B2) results in a second independent equation 
on the first order quantities. Therefore, the first-order 
inner equations for the perturbation quantities can be taken 
to be the set consisting of the first of Eqs. (Bl) and (B2) 

2 

and the difference between the two 0(e ) equatxons: 


■rix = “ ■ 


X X 
- U-,r, 


(B3) 
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APPENDIX B 


Finally, the first equation is used to eliminate from 

Eq. (B3) to give the first-order inner equations in the form 


n 


(B4) 


1 1 . ro 1 

Ult + [2y^ 


(Y-1) r^ - 2m^ ] 


2aX i 

I^^IX (Y+Dm^ ^^^1 


(Y-l)r^] = 0 


Equations (B4) are the lowest order nonlinear equations of 
motion which govern sound propagation in the vicinity of a 
near-sonic throat. 

A second order equation can be found by forming the second 
order Riemann invariant (Eqs. (3.11)) 


P2 = P2 + 


r^ - 1 (r^)^ - m 
^2 4 ^1^ 1^1 


(B5) 


Using the second of Eqs. (Bl) and the relation ” 0' 

which is implied by Eq. (B4) yields 



3 

3X 




Y-3 / i. 2 „ i> 


0 


and hence = P 2 (t) . The relations = 0, P 2 = ^2^^^ 

the second of Eqs. (B4) are equivalent to the equations derived 
in Section IIIB. In fact the second of Eqs. (B4) can be trans- 
formed to Eq. (3.14) by letting = y^-r^. 

The solution for y^ and r^ can be expressed in terms of tt 

by use of the relations Q, = 2yw y^ Eqs. (3.15) . 

This yields 

y^ = - [7T(X,t) + m^(X)] . (B6) 

Using Eq . (3.20) for it gives the leading order inner solution 

y^ = - rj = :^ [-(aX^ + c|)^(t))^/^ + (aX^ + 1)^^^] , (B7) 
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Where T(X,t) is given implicitly by the characteristic equation 
( 3 . 21 ) . 
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Appendix C. FOURIER EXPANSION OF THE INNER SOLUTION 

In Section IIIC the Fourier series for the outer expan- 
sion of the inner solution was found assuming the Fourier 

2 

expansion of the function 0 (x) was known. In this appendix 

2 

the details of finding this expansion of cj) [x(x,t)] will be 

2 

provided. Since tf) is a period function of period 2iT/fi in t 
a series in the form 


4 > 


2 



+ 


1 


n=l 


+ + 

C (x) cos nJ^t + D (x) sin nUt 
n n 


(Cl) 


will be sought, with 


f + 1 

^n 

< 

-f- 

q 

IT 

■2it/Q 

(p^ (t (x,t) ) • 

cos f^t 

> 

D" 

. 


sin Qt 

n 

V. J 

u 



(C2) 


+ + 

It should be noted that C and D are functions of the outer 

n n 

variable x since the outer limit of the inner solution is 
being considered. The expression (C2) for the coefficient 
will now be simplified and written in the form suitable for 
matching by explicitly exhibiting the x dependence implicit 
in Eq. (C2) . Differentiating the characteristic relation 
(Eq. (3.29)) 


t 


T 


+ 


sgn X 

1/2 

a 


In 


1 (i) 


+ g(x;e) 


(C3) 


yields 


dt 

dx 


1 + 


sgn X <p ' 


J72 


4 > 


(C4) 


The variable of integration in Eq. (C2) is changed from t to 
X by using Eqs. (C3) and (C4) yielding 
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D' 


n 

+ 

)~ 

n 




TT 




X < 


COS nfi 

sin nfl 


|j + g(x; e) + -^2^7^ In |<()(t)|J 


mr 


Here the fact that 4>(t) is periodic with period 2n/n has 
been used to determine the limits of integration. Integrating 
the (() ' term in C~ and D“ by parts yields 


_± nJ2 _ + 

C_ T-Vo sgn X D 


" 2aT^ 


and 


D ■! 1 

n 0-.1/2 


n 


TT 


2n/Q 


(p^(r) cos nfl |j + g(x;e) + ^^^ 2 ^ I ij (C 5 ) 


sgn X C 


2a 


n 


U 

IT 


2tt/J^ 


(j)^(T) sin ni^ T + g(x;e) I 4* I dr 


(C 6 ) 


0 


respectively. Multiplying Eg. (C6) by i and adding it to 
Eg. (C5) yields 


tW n D I + i 


C^j = ^axp (ln!Jg)I* (C 7 ) 


+ + 
where I was defined in Eg. (3.32). Since I is independent 
n ± ± ^ 

of X the X dependence of and has been explicitly found 


in the right side of Eg. (C7) . 
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In order to solve Eqs. (C6) and (C7) for the Fourier 

a. a. I 

coefficients 


C and D the term 
n n 

1“ 

n 

is written as 



+ 

+ + , + 1 
I = E + iH" = 

I" 1 

i tan ^n 

n n n ' 

n ' 

^ + 
TT 



n 


Then after a straightforward calculation the expressions 

(3.31) are found for C~ and D~. 

n n 
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Appendix D. THE CROCCO-TSIEN SOLUTION 

An exact analytical solution of the outer acoustic 
equations (2.25) was discovered by Tsien and applied extensively 
by Crocco to problems of combustion instability in rocket engines 
(Ref. 10). Their solution was for a basic steady flow that 
was subsonic upstream of the throat and supersonic downstream 
of the throat. This solution was generalized in Ref. 6 to 
totally subsonic steady flow. This generalized solution is 
presented here. 

Suppose the Mach number distribution (x) of Eq. (2.17) 
is of the form 

- 1/2 

Mq = (1 - k|x| ) (1-k1x|)^] |x| (Dl) 

where K is a constant. Expansion of Eq. (Dl) for small x 
yields 

Mq (x) = 1 - |x| + I (Y^ - DK^x^ + . . . . (D2) 


This is expansion of the form of Eq. (2.17) with 

2 ^ 2a 
^ (Y+1) 

and 

^20 = f 


The area ratio a (x) corresponding to Eq. (Dl) can be deter- 
mined from Eq. (2.13); it is shown in Figure 4 for three 
values of the parameter K over a range 0 ^ x ^ 0.5. Tsien 
and Crocco showed, by introducing a new independent variable 

z = (1-k|x1) into Eqs. (2.25) that these equations were 
equivalent to 


,2 . 

,, , ^On ^ . ig T _ On 

z (1-z) 7 , 2 [1 + z 


dz 


Y+1' 


dz 


, i3 (2+2i) . 

^ 2 (Y+1) On 


0 


(D4) 


dr 


ton = [(Y-l+i3)ro^ 


- (Y+1) (1-z) 


dz 


[2+iB]"^ 


(D5) 
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where 


B = -(l±i) 


1/2 


K 


sgn X 


(D6) 


Eq. (D4) a hypergeometric equation with complex coefficients 
Ref. 20 and thus its general solution can be written as 


l-Si 


^In^^ f (“5^/~j/2 £;1 z) + E 2 j^E ( j ,k, 5.; 1 — z) , 


(D7) 


where 


^ - ^-1 1 a 

D 2 ~ 2 ^ ’ 


k = »^ + 1 A , 


1 = 2 + 


2i3 

Y+1 


A = [1 + ^ 3^] 

(Y+1) 


1/2 


and F is the standard hypergeometric function of the indicated 
complex arguments. Calculating r^ from Eq. (D7) by using the 

well known properties of the hypergeometric function gives 

from Eq. (D5) in the form 


‘On 2+i3 \ 2n 


e; 


[ (y-l+i3)F^ + (1 -z)F2] 


1-5., 


+ E£^ [(2Y+i3-(Y+l)^) (1-z) F 3 


(D8) 




where 


F^ = F( j ,k,5,;l-z) , F^ = F ( j +1 ; k+1 , £+1 ; 1-z ) 


F = F(-k,-j,2-5,;l-z) , F, = F (-k+1 , - j +1 , 3-5, ; 1-z ) 


Finally, if r^^ in Eq. (D7) is multiplied by exp infit and 
expanded for small values x (i.e., z-1 -> 0) it must be of 
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the form given in Eq. (2.45)/ which determines E^ and E» in 

terms of A7 and B7 : 

In In 



(2K)^ ^ 



and 




(D9) 


If the anechoic boundary condition 


P^^(-x ,t) = 5 
e 

is applied then using Eqs 
Eq. (3.47) namely 


y (M_y^ + r 

0 On 

n=0 

(Dl), (D7), 


)^inat . 

On 

(D8) and (D9) yields 


B, = (2K)^ ^(r (-X ))"‘^a7 

In n e In 


where 


n e 


(2+iB)"^[(Y-l+iB)F^ + 


Mo(-Xe) 


(2+ie ^) 


(2y+i3-(Y+l)S-) (l-z^)^~S3 + ^^ 4 ] 




(DIO) 


and 


^e = 


K -x^|] 


In this expression the hypergeometric functions are evaluated 

= z . 
e 

If the anechoic boundary condition 


at z = z . 

e 


q"(x t) = « I «o“on- = “ 

n=0 
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is applied then a similar procedure to that given above yields 

where is given by an expression similar to Eq. (DIO), 

with the following changes: the signs of terms divided by 

1/Mq should be changed and all quantities in (DIO) should be 

evaluated at x = x > 0. 

0 
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FIG. 13. Percentage dist: 
total power amoi 
for varying Mac] 
= 1, a = 1, N 
= .03. 
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FIG. 16. Percentage distribution of total power 
among harmonics for varying throat 
curvature; = .01, n= 1, e = .1. 
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FIG. 18. Percentage distribution of total power among harmonics 
for varying source strength (downstream power flow) ; 
i^ = l, a = l, e = .1. 
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G. 20. Percentage distribution of 
for varying throat curvatui 
N = .01, Q. = 1, e = .1. 
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